Quantitative description of Josephson-like tunneling in z/^ = 1 quantum Hall bilayers 
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At total filling factor vt = 1, interlayer phase coherence in quantum Hall bilayers can result in a 
tunneling anomaly resembling the Josephson effect in the presence of strong fluctuations. The most 
robust experimental signature of this effect is a strong enhancement of the tunneling conductance 
at small voltages. The height and width of the conductance peak depend strongly on the area and 
tunneling amplitude of the samples, applied parallel magnetic field and temperature. We find that 
the tunneling experiments are in quantitative agreement with a theory which treats fiuctuations 
due to meron excitations phenomenologically and takes tunneling into account perturbatively. We 
also discuss the qualitative changes caused by larger tunneling amplitudes, and provide a possible 
explanation for recently observed critical currents in counterflow geometry. 



I. INTRODUCTION 



The existence of a bilayer quantum Hall state at 
total filling factor = 1 is well-established both 
theoreticalljii"— and experimentallji^i"— . This state is 
characterized by remarkable electronic properties such as 
counterflow superconductivity and a Josephson-like en- 
hancement of tunneling between the two layers>ii^"— Its 
formation is controlled by the relative magnitude of intra- 
and interlayer Coulomb interactions, and therefore de- 
pends critically on the ratio of interlayer separation d and 
magnetic length Ib = \/h/eB. At large d/ls, there is no 
quantum Hall effect and the bilayer system behaves quali- 
tatively like two independent composite fermion systems. 
On the other hand, at small d/ls, the interlayer Coulomb 
interaction induces an exotic quantum Hall state, which 
can be described as an exciton Bose-Einstein condensate 
or as a pseudospin ferromagnet.— i^i^ Here, the pseudospin 
is formed from the two- valued which layer quantum de- 
gree of freedom, and condensation occurs because the 
electrons can lower their interlayer exchange energy by 
entering a state with uncertain layer indexii 

Experimentally, the most spectacular effect arising due 
to the interlayer phase coherence at small d/lg is a huge 
enhancement of the tunneling conductance at small volt- 
ages. In the incoherent state at large d/lB, the tun- 
neling in quantum Hall bilayers is strongly suppressed 
by a Coulomb gap,—"— which is caused by strong elec- 
tronic correlations and the squeezing of wave functions 
by the magnetic field rather than by small geometric 
dimensions 1^ Namely, in a strong magnetic field each 
layer is in a correlated state of its own, and the tunneling 
event can be considered as an electron suddenly brought 
from one layer to the other, resulting in an excited state 
of the system. Relaxation towards the ground state takes 
place when the energy is carried away by collective ex- 
citations, and this excitation energy must be provided 
by a sufficiently large external voltagei^ On the other 
hand, if the layers are closer together, the Coulomb gap 
decreases^ due to interlayer correlations^, and finally at 
some critical value of d/ls a strong and sharp peak ap- 



pears in the differential conductanc o^^i'^^ , indicating the 
transition to an interlayer coherent state with a quantized 
Hall drag and small counterflow resistances^Sii^i^"— 
The coherent state was further characterized by deter- 
mining the dispersion relation of the collective Goldstone 
modo^^ and the degree of spin polarizatio n'^^i^^i^^ . More- 
over, the phase-transition from incoherent to coherent 
stat o^^i^-^i'^^i^^i'^-^i^^ was studied in some detail. 

The Josephson effect with a zero bias supercurrent can 
be observed when there is a phase difference between 
two superconductors separated by a tunnel junction. In 
a quantum Hall bilayer system, the individual layers are 
not superconducting by themselves, and only the two lay- 
ers combined can exhibit phase coherence. However, in 
the coherent state of a clean system with tunnel coupling 
between the layers, a deviation of the pseudospin orien- 
tation from the minimum energy direction is expected 
to give rise to a zero bias tunnel current between the 
layers^, in analogy to the Josephson current in a super- 
conducting tunnel junction. Experimentally, instead of 
a zero-bias supercurrent, a Josephson-like enhancement 
of tunneling^^i^i^"— i^"— at small interlayer bias volt- 
age was observed. A significant amount of theoretical ef- 
fort has been devoted to understanding the finite height 
and width of the conductance peak of this Josephson- 
like tunneling at small voltagesiSrJ^iiiii^iS^ as well as the 
magnitude of the critical tunneling curren t^^i^^i^^ . These 
approaches differ in how the bias voltage, quasiparticles 
and disorder are introduced into the theory. In one type 
of approacbi2i2ii^, the clean limit without charge disor- 
der is considered. Then, the critical current is given by 
the maximum current at which the order parameter can 
be static, and the tunneling conductance cannot be de- 
scribed in terms of condensate dynamics alone, because it 
also depends on the fermionic quasiparticles which trans- 
port the charge between leads and bulkJ^i^^ 

Experimentally investigated samples are probably not 
in the clean limit, and experimental results indicate that 
disorder-induced topological defects, so-called merons, 
are important .^iS^—'i^^-i^'i^'^'S^'^'S^ Merons carry a 
charge ±e/2 and are characterized by their vorticity 
and the layer in which the charge resides J^iiii^ Be- 
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cause merons are nucleated by charge disorder, they ex- 
ist also at low temperatures, and their dynamics is ex- 
pected to give rise to dissipation and to strongly affect 
the interlayer tunnelingj^ ^^^i^^^^^'^^'^^i^^i^"^ Merons can 
be taken into account in the description of Josephson- 
like tunneling by introducing a phenomenological vortex 
fieldJ^-i2>^>25 Boi-Ij dynamioiaJi and statioH>24^ cor- 
tex fields have been considered in the literature. Despite 
similar starting points, static and dynamic approaches 
predict different characteristic features for the tunneling 
current. 

Spatial fluctuations in the vortex field are accounted 
for by a correlation length ^, which governs the decay of 
the pseudospin ferromagnetic order. Spatial fluctuations 
are caused by the formation of compressible puddles, 
which arc thought to arise due to local fluctuations in the 
density of dopants. If the local density of merons created 
by charge disorder is sufficiently high, the merons can 
screen the random potential and delocalize. As merons 
are vortices of the order parameter field, spatial corre- 
lations in the order will decay on the typical distance 
between puddles, which in turn is determined by the set- 
back distance to the dopant layer, such that one expects 
^ ~ 100 - 200 nmlSi^. 

In the present manuscript, we study the influence of a 
dynamical vortex field on interlayer tunneling deep inside 
the coherent phase. Dynamical fluctuations of the order 
parameter are caused by the dynamics of merons in this 
picture. At temperatures lower than the energy gap but 
still comparable to it, one can imagine that thermally 
activated hopping of merons from one puddle to another 
is the dominant source of fluctuations. By mapping the 
quantum Hall bilayer to a classical two-dimensional XY 
model with a symmetry-breaking field and with disorder, 
Fertig and Straleyi^ have found that disorder nucleates 
strings of overturned spins, which connect vortices and 
antivorticcs at their ends. At low temperatures, this state 
shows glassy features and gives rise to anomalously large 
fluctuations of the vortex field. The temperature depen- 
dence of interlayer tunneling in a coherence network of 
puddles and ordered regions was analyzed in Ref. fisl . 

We model dynamical vortex field fluctuations by in- 
troducing an exponential time decay of local pseudospin 
corrclationsiSiii, which is governed by a correlation time 
T^. In this approach, tunneling between the two layers 
can be treated perturbatively, and Josephson oscillations 
in the presence of vortex field fiuctuations give rise to the 
finite tunneling peak observed in the experiments. This 
effect strongly resembles the Josephson tunneling in small 
Josephson junctions in the presence of thermal or quan- 
tum fluctuations^"—. We find that the current-voltage 
(I-V) characteristic of interlayer tunneling is character- 
ized by a voltage scale 



Vq = hu/e^ 

and by the scale for the maximum tunnel current 
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Here, u is the velocity of the pseudospin wave mode, 
'-^SAS is the tunnel coupling between the two layers, is 
the area of the sample and Ps is the pseudospin stiffness. 
The functional form of the I-V characteristic is controlled 
by the temperature dependent decoherence rate 



a(T) 



UT^{T) 



(3) 



which for instance enters into the zero bias conductance. 
For small values of a we obtain 



Go(r) 



/q 2 
a{T) 



(4) 



While the zero bias conductance increases when a — ?> 0, 
the maximum tunnel current saturates in this limit to 
a value tt/q. An in-plane magnetic field suppresses the 
conductance peak, and we find that the characteristic 
field scale for this suppression is ^o/{d£^), where $0 ~ h/e 
is the fiux quantum. 

We compare our theoretical predictions to the experi- 
mentally observed dependencies of the tunneling current 
on sample area, tunneling amplitude, applied parallel 
magnetic field and temperature. We find that experi- 
mental results and theoretical predictions are in excel- 
lent agreement. Despite the comparatively large number 
of parameters in our theory, we find that all parameters 
can be determined uniquely from experimental data and 
that the parameter values obtained in this way are in very 
good agreement with theoretical estimates. We find that 
the apparent temperature dependence of the tunneling 
conductance^ Go oc exp(— T/Tq) can be explained as a 
crossover between an activated and a power law temper- 
ature dependence of a(T), which in turn has a natural 
explanation in terms of vortex field fiuctuations. 

More recently^"—, intriguing observations were made 
in bilayer samples with larger tunneling amplitude. In 
particular, jumps of the tunnel current were observed and 
described in terms of critical tunnel currents about two 
orders of magnitude larger than in previous experiments. 
We find that these observations can still be explained in 
the framework of a perturbative approacbiSili when tak- 
ing into account that the tunneling resistance between 
the layers is so small in these samples that the other 
resistances in the system, such as circuit, contact, lon- 
gitudinal and Hall resistances, become important. Our 
approach also provides an explanation for recently ob- 
served critical currents in a counterflow geometry^. This 
surprisingly good agreement with the experiments over a 
wide range of tunneling amplitudes raises questions re- 
garding the region of validity of the perturbative treat- 
ment of the tunneling. We give a short discussion of this 
important theoretical issue in Section IVTl 

Our paper is organized as follows. In Section |ll] we de- 
scribe the theory of interlayer tunneling in the presence of 
fiuctuations caused by the meron excitations. In Section 
lllll we compare the experimental tunneling I-V character- 
istic obtained using samples with extremely small tunnel- 
ing amplitude to our theoretical predictions. We find a 
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good quantitative agreement between theory and exper- 
iments. In Section IIVI we apply the theory to describe 
the quahtative changes observed in recent experiments 
where the tunnehng amphtude of the samples was signif- 
icantly larger. In Section |Vl we compare with previous 
theoretical approaches, and finally in Section [VII we give 
a summary of our findings and discuss the limits of the 
validity of our approach. 



II. THEORY OF INTERLAYER TUNNELING IN 
THE PRESENCE OF MERON EXCITATIONS 

We concentrate on the tunneling deep inside the coher- 
ent phase at vt = noh/{eB) = 1, where no = {2iTlg)^^ 
is the average density. We assume that the real spin is 
fully polarized in agreement with the experiments^'^'^. 
The pseudospin is therefore the most relevant degree of 
freedom in our problem: The low-energy excitations of 
the system are the pseudospin waves and the merons. 
The theoretical analysis presented in this section largely 
follows that of Stern et al. in Ref. [l^. In the absence 
of tunneling, the pseudospin waves are described by a 
Hamiltonian densityii^iii^ 



2r 



(5) 



where to = (cos ip, sin tp, ruz) is the pseudospin vector, ps 
is the pseudospin stiffness and T is the capacitance per 
area. The momentum conjugate to (/? is tt = hnQmz/2. 
Physically the first term in equation ([5]) arises from the 
loss of optimal Coulomb exchange energy: If the pseu- 
dospin varies in space, the spatial part of the many- 
particle wave function cannot be fully antisymmetric. 
The second term measures the capacitive energy. Be- 
cause of exchange effects between the layers F is expected 
to be strongly enhanced from its electrostatic valueJ^^^^ 
The dispersion relation for collective pseudospin waves 
(Goldstone mode) is 



ujj: = uk, 



where the velocity of the pseudospin waves u is 




(6) 



(7) 



The tunneling energy should be included in the Hamil- 
tonian ^ as T+-|-T_— , where the tunneling operators 



^±ieV t / h ^±iQ B X 



(8) 



Here A = l^sAS /{^t^I'b)^ Qb = eB^d/H, is the in- 
plane magnetic field and V is the interlayer voltage, 
which we assume to be constant in space. The merons 
are included phenomenologically with the help of the vor- 
tex field ji2r2^i2i Using the commutation relation for 



ip and TT it is easy to verify that the tunneling operators 
T± change the charge on the capacitor plates by ±e. 

Both counterflow and tunneling experiments suggest 
that at least part of the merons are mobile. Therefore, we 
assume that the vortex field is dynamic. The tunneling 
current depends on the details of vortex field dynamics, 
which we assume to be characterized by a correlation 
length ^ and a correlation time t^. For the correlation 
function we make the specific assumption 



-r/ip-\t\/T^ 



(9) 



which appears to produce a good agreement between the- 
ory and experimental results. By studying the experi- 
mentally found dependence of the tunneling current on 
a parallel magnetic field, we find strong evidence that 
^ ^ 100 — 200 nm, in agreement with earlier theoretical 
estimate o^°'^° . Because pseudospin fluctuations are un- 
correlated on scales much larger than ^, we can at least 
qualitatively think of the small domains of size as in- 
dividual Josephson junctions. These small domains are 
well-coupled to each other by counterflow and transport 
currents. We assume that the tunneling can be consid- 
ered perturbatively due to the reasonably small corre- 
lation time caused by the strong fluctuations of the 
vortex fieldi^. The validity of the perturbative treat- 
ment of the tunneling has been well-established^ in small 
Josephson junctions in the presence of fluctuations, but 
the justification of similar assumption in quantum Hall 
bilayers is questionable especially at low voltage o^^i^^ es- 
sentially because the phase (p describes both the counter- 
fiow currents and the tunneling simultaneously. Never- 
theless, we find that the perturbative treatment of the 
tunneling gives a good quantitative description of the 
Josephson-like tunneling for all voltages and a surpris- 
ingly large range of tunneling amplitudes. We discuss 
this unexpected result and the limits of validity of our 
approach in more detail in Section IVTl In this section we 
concentrate on a situation where the tunneling resistance 
is larger than the intralayer resistances, Rxx and Rxy, so 
that the interlayer voltage can be considered homoge- 
neous. However, we stress that our approach can easily 
be extended to a case where the interlayer voltage varies 
slowly in space. We believe that the interlayer voltage 
is very homogeneous in the experiments where the tun- 
neling amplitude is small^^i^i^"— i^^"— . In the case of 
larger tunneling amplitudes the inhomogeneities will af- 
fect the details of the tunneling I-V characteristics, but 
we believe that they still do not play an important role in 
the main qualitative features observed in the experiments 
reported in Refs. list - liel (see Section |IV|. 

The I-V characteristics can now be calculated using 
Eqs. dS])-® and Fermi's golden rule, similar to the case 
of Josephson junctions^i^. The result of this calculation 
is 

I{V,B.) = 1^ Hdi I d^r G™(r,t)e-^(-*)/2 



Jo J 

■ C{r,t) . eVt 

X sm — - — cos\Qbx) sm — — , 
^ n 



(10) 
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where 

D[r, t) = j2 — ^ cos(k • r) cos[ukt)\ cotn — ^ — , 



hu 



sin{ukt) cos{k ■ r) 



2-Kp, 



-B{ut - r) 



2 1-1/2 



(11) 



(12) 



/3 l/fcijT and G™(r,t) = (e*'^™(r,t)g-»<p„(o,o)^ ^j^g 
correlation function for ipm- This result has been found 
previously in Ref. [lol . but for completeness we give a 
detailed derivation of equations ([T0| - (fT2|) in Appendix 

m 

In order to calculate we need to define an "ultra- 
violet cut-off momentum"— /cq = k\/2/^_b, where it is 
expected that k « 1. Although some of the parameters, 
like ps, r and AgAS, can be theoretically estimated^iiSi^i, 
there remain reasonably large uncertainties in their ac- 
tual values. We find that it is possible to determine most 
of the parameters based on experimental data and the 
obtained values arc in good agreement with the theo- 
retical predictions. Before comparing theory with ex- 
periments, we will make certain simplifying assumptions. 
The theoretically estimated value for pseudospin stiffness 
is ps ~ 0.4 K— , and according to Ref. [2^ this parame- 
ter is likely to be the most reliably known continuum 
model parameter. Using this value of Ps and the values 
of the other parameters found below, we have numerically 
checked that D depends only weakly on temperature in 
the experimentally relevant temperature range. There- 
fore, we can use everywhere the zero temperature value 



huka 
2Trps 



2Tr^p,Th 



(13) 



For _B = 2 T, we obtain Dq '--^ 3.3, but an accurate es- 
timate is not possible due to the uncertainty related to 
the ultraviolet cut-off momentum. With similar assump- 
tions, we have numerically checked that we can use the 
following first order expansion in Eq. (jlOp 



sm ■ 



C(r,t) C(r,t) 



(14) 
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These simplifications were justified also in Ref. [l^. 

By using the Eqs. ([T0| -(|9|) and the approximations out 
lined above, we get {R ~ r/^ and q ~ k^) 



Iq / dq 



+ {V/Vo - qY a2 + [V/Vo + qf 



j dRR e~^ Jo{qR)MQBiR), 



where 



'Do/2 



h Ps 



/4 



647r2 



(15) 



(16) 



and Vq and a are given by Eqs. ([T]) and ([3]), respec- 
tively. In Appendix|B]we give an alternative derivation of 
Eq. ()15p where the quantum fluctuations are neglected. 
Similarly as in Ref. Il2l . we solve the Sine-Gordon equa- 
tion perturbatively, but in contrast to Ref. IT^ we also take 
into account the dynamics of the vortex field. Finally, we 
obtain the same expression ([T5|) for the TV characteris- 
tics except that now the factor e"^"/^ docs not appear 
in the equation for Iq. We interpret this factor as the 
effect caused by quantum fluctuations. 

In the absence of a magnetic field, = 0, we obtain 



I{V) = h dq 



+ iyiv^ - qf 



O? + {VIVq + qY 



(1 + g2)3/2 • 



(17) 



In the limit V 
conductance Gq 



and for small values of a, the peak 

= #ly=o is gi^^"^ by ©■ The 
maximum tunneling current in the limit a — > saturates 

to a value tt/q. 

We expect that Vb, /o and ^ depend only weakly on 
temperature. The strong temperature dependence ob- 
served in experiments would then be caused by the tem- 
perature dependence of T^p. In our opinion, this assump- 
tion is in agreement with more detailed theoretical mod- 
els where the meron dynamics were analyze d^^'^^'^^ . We 
expect that for temperatures comparable to the thermal 
activation gap Ag, the temperature dependence of 
arises mainly due to thermally activated hopping of the 
merons. On the other hand, the temperature dependence 
of at temperatures much smaller than Ag could origi- 
nate from glassy excitations of the pinned meron system, 
which were predictedi^ to act as a dissipative environ- 
ment. 

The determination of the different theoretical param- 
eters of the model is discussed in detail in Section IIIII 
We summarize the methods used to determine them and 
the values we obtain in table HI The values of the param- 
eters obtained using the experimental data are in good 
agreement with the theoretical expectationsiiiSi2£i2i. 



III. QUANTITATIVE DESCRIPTION OF 
JOSEPHSON-LIKE TUNNELING IN THE CASE 
OF SMALL AsAS 

We start by considering the first experimental observa- 
tions of the Josephson-likc tunnelin g^^'^^i'^^'^^ . Our the- 
ory should be most applicable to these experiments, be- 
cause the tunneling amplitudes of the samples were ex- 
tremely small, and therefore higher order tunneling pro- 
cesses can be safely neglected. In addition, the tunneling 
resistance is larger than the intralayer resistances Rxx 
and Rxy , so that the interlayer voltage can be considered 
to be constant in space. We expect that in samples with 
sufficiently small As AS tunneling is a bulk phenomena 
which takes place homogeneously throughout the whole 
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Parameter 


Method of determination 


Value 


Ps 


Theoretical estimate 


~ 0.4 K 




Dependence of the tunnehng 
conductance on B\\ 


~ 130 nm 


u 


Resonant enhancement of the 
current at eV = hnQs 


^ 14 km/s 


r 


Using estimates for ps and u 


^ lOeeo/d 




Height and width of the 


Temperature 




conductance peak 


dependent 


AsAS 


Tunnehng I-V characteristics 


~ 10 ^JK 




at B = 


~ 100 mK 



TABLE I: Determination of the parameters. The different 
estimates for Asas ~ 10 and As as ^ 100 /iK correspond 
to samples considered in Sections IIIII and II VI respectively. 



sample area, independent of the sample geometry. Re- 
cently, it was found^ that the peak in the tunneling 
conductance is proportional to the area of the sample. 
Moreover, it seems that in these experiments the tun- 
neling current at all voltages scales proportionally to the 
area. These observations are explained naturally within 
the present approach, which predicts the current scale /q 
to be proportional to the area [see Eq. ([H])]. 

The layer separation and the surface area in these 
experiments2^i^i^i24 are d = 27.9 nm and = 250 • 250 
/im^, respectively. The tunneling amplitude can be 
determined^ by solving the Schrodinger equation for the 
specific structure. Another possibility is to use the rea- 
sonably well-understood^2i^ I-V characteristics at i? = 
as explained in Ref. |^. These two methods are in good 
agreement with each other—, and therefore we are rea- 
sonably confident that Ag^s in these experiments is ap- 
proximately 5 — 10 ^K. We discuss the experimental and 
theoretical estimates for the other parameters ^, and 
u below. 



A. Dependence on parallel magnetic field 

The in-plane magnetic field in quantum Hall bilayers 
has two main effects on the Josephson-like tunneling. It 
suppresses the small-bias current and results in resonant 
enhancements of the tunneling current at voltages satis- 
fying eV = HuQb (see Fig. [T]). In this section we dis- 
cuss the physics behind these effects and determine the 
values of ^ and u by comparing the theory and exper- 
iments. We concentrate on the low voltage part of the 
I-V characteristics where the experimentally measured 
[Fig. [1] (a)] and theoretically calculated [Fig. [T] (b)] tun- 
neling currents are in good agreement with each other. 
The theoretically calculated tunneling currents at large 
voltages are typically significantly smaller than the ex- 
perimentally measured currents [cf. Figs. [1] (a) and (b)]. 
The origin of this discrepancy is probably an incoher- 
ent contribution to tunneling discussed in more detail in 
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FIG. 1: (a) Measured I-V characteristics for d/lB = 1.61, T = 
25 mK and different values of parallel magnetic field _B|| =: 
0,0.11,0.15,0.20,0.23,0.29,0.35,0.41,0.47 T. (b) Theoretical 
I-V characteristics obtained using Eq. (|15[) for Vb = 73 /^V, 
la = 7.1 pA, a = 0.04, ^ = 126 nm and the same values of 
parallel magnetic field. The determination of the theoretical 
parameters is described in the text. The experimental curves 
are reproduced from the experimental data provided by I. B. 
Spielman et al. and reported in Ref. [3^ . 



Section ED 

As already pointed out the domains of size ^ can, 
at least qualitatively, be considered as individual Joseph- 
son junctions, which are well-coupled to each other by 
counterflow and transport currents. If a magnetic flux 
<& is applied across a Josephson junction, it tries to in- 
duce a circulating current inside the junction)^ Due to 
this effect, the critical tunneling current of the Josephson 
junction oscillates as a function of $/$o, where $o is the 
magnetic flux quantum. The current completely vanishes 
if an integer number of flux quanta is applied across the 
junction and the dependence, in general, is typically re- 
ferred to as Fraunhofer diffraction pattern by the analogy 
with the case of light passing through a narrow rectan- 
gular sliti^ Based on our analogy, we can now define for 
each domain a magnetic flux penetrating the area be- 
tween the layers as $ = B\\£^d. Similarly to the case of 
Josephson junctions, we expect that the current vanishes 
quickly as a function of $/<I>o- The important difference 



B,=0.47T 
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FIG. 2: Experimentally measured and theoretically calculated 
small bias tunneling conductance as a function of parallel 
magnetic field _B| | . The experimental results were obtained 
at r = 25 mK. The small-bias tunneling conductance was 
calculated using Eq. (fT5|) for a = 0.04 and ^ — 126 nm. The 
value of ^ was obtained by fitting Lorentzians to the theoreti- 
cal and experimental results and by setting the widths of the 
Lorentzians equal to each other. Inset shows the theoretical 
small bias conductance as a function of Qs^ (dots) and the 
corresponding Lorentzian fit (line). The experimental data 
was provided by I. B. Spielman et al. and has been reported 
in Ref. [H. 



to standard Josephson junctions is that in our case, the 
size of the junctions is not wcU-dcfincd and therefore the 
Fraunhofer diffraction pattern is "washed away". Here 
the exact form of the r dependence of the vortex-field 
correlation function Gm{r,t) also becomes important. It 
defines the size distribution of our Josephson junctions 
and therefore it strongly affects the shape of the small- 
bias tunneling conductance plotted as a function of i3| | . 
We find that the exponential dependence given by Eq. ^ 
gives significantly better agreement with the experiments 
than the Gaussian dependence suggested in Ref. 

Figure [5] shows the experimentally measured small bias 
tunneling conductance as a function of i?| | and the corre- 
sponding theoretical result obtained using Eq. Both 
dependencies are approximately Lorentzian as demon- 
strated for the theoretical curve in the inset of Fig. O The 
half width at half maximum (HWHM) of the Lorentzians 
fitted to the experimental and theoretical curves are 
^HWHM ^ 0.138 and QgWHM^ ^ o.735. Using the rela- 
tion Qb ^ eB\\d/h, we find ^ w 126 nm. This method 
of determining ^ is extremely robust against variation of 
the other theoretical parameters. The obtained value is 
completely independent of the values of /q and Vg , and it 
changes only by a few percent if the parameter a is varied 
between 0.001 and 0.1. Moreover, the value of ^ obtained 
in this way is in excellent agreement with the theoreti- 
cal expectations in Refs. IToll20ll23 . It is also important 
to note that in clean samples a Fraunhofer pattern has 
been predicted^, and much smaller magnetic fields will 
induce circulating tunneling currents. Therefore, the in- 
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FIG. 3: (color online), (a) Experimental measurement of 
(fl/dV^ for different values of parallel magnetic field -B|| —: 
0.20,0.23,0.29,0.35,0.41,0.47 T. (b) Theoretical curves ob- 
tained using Vo = 73 fiY, lo = 7.1 pA, a = 0.04, ^ = 126 nm 
and same values of B\\. Dots show the expected locations of 
the resonances V* = huQB/e, which for sufficiently large Qb 
are very close to the local extrema of d?I/dV^ . In theoreti- 
cal curves we have introduced larger offset (2 ■ 10"'^ Q.~^V~^) 
between the curves than in experimental results (0.8 ■ 10""^ 
Q,~^V~^). The experimental curves are reproduced from the 
experimental data provided by I. B. Spielman et al. and re- 
ported in Ref. [s^ . 



plane magnetic field dependence of the small-bias tun- 
neling conductance clearly manifests the important role 
of the merons in the Josephson-like tunneling. 

The second interesting effect caused by the in-plane 
magnetic field are the additional enhancements of the 
tunneling current at voltages V*{B\\) satisfying eV* = 
HuQb. The origin of these resonances is nicely explained 
in Refs. ITol - [T2l . The parallel magnetic field allows only 
tunneling between states that differ by a momentum Q b , 
and energy conservation requires that the energy of these 
states differs by eV . In the absence of a dissipative en- 
vironment, there would be just a single linearly dispers- 
ing collective mode = uk and therefore both condi- 
tions can be satisfied only if the voltage satisfies the res- 
onance condition eV* = HuQb. Therefore the resonant 
enhancements of the tunnehng current in the I-V char- 
acteristic essentially map the dispersion relation of the 
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pseudospin wavesJ^— Because of the finite correlation 
length and time caused by the merons, these resonances 
are not sharp but they can still be clearly seen both in 
the experimental and theoretical curves by plotting the 
second derivative (Pl/dV^ as a function of V as shown in 
Fig. [3l The second derivative is not sensitive to any lin- 
ear background conductance and therefore it shows more 
clearly the nonlinear features in the I-V characteristics. 
In particular the voltages V* at the resonances can be 
determined by finding the extrema of (Pl/dV^, which 
indicate a maximum of the curvature of the I-V charac- 
teristic. This result was already anticipated in Refs. [2^ 
andlH, where the velocity u of the pseudospin waves was 
obtained experimentally by determining for each mag- 
netic field the voltage V* where (Pl/dV^ has an ex- 
trcmum, and by fitting a linear relation V* = huQB/e 
to the data points. Wc have theoretically confirmed the 
validity of this approach by numerically determining the 
locations V* of extrema of d^I/dV^ for different values 
oi Qb- For Qb^, ^ 2, the numerically determined func- 
tion V*{Qb) can be approximated by the linear relation 
V* = huQs/e within an accuracy of few percents. In 
addition, the numerically determined function V*{Qb) 
is completely independent of Iq and practically indepen- 
dent of a, i.e. it does not change within our numerical 
accuracy when a is varied between 0.001 and 0.1. There- 
fore, the experimentally found velocity u « 14 km/s can 
be considered as a reliable estimate of the actual veloc- 
ity of pseudospin waves. This value of u is somewhat 
smaller than theoretically expected^'i^'^'^^i^, but prob- 
ably still within the theoretical uncertainties because the 
parameters ps and F are not known accurately. 



By using our estimates for the correlation length ^ and 
Goldstonc mode velocity u, we find Vb ~ 73 ^V, which 
yields good agreement between theoretical and experi- 
mental curves in Figs. [T] and [3l We can now determine 
the parameters a and Iq by using the measured values of 
small-bias tunneling conductance Go — 4.68 • lO^*' il^^ 
and maximum tunneling current /max — 18 pA at _B|| = 
0. We obtain a « 0.04 and Iq « 7.1 pA. By assum- 
ing that Ps « 0.4 K, we estimate using Eq. ([T5| that 
Iq ^4—16 pA (corresponding to Asas = 5—10 pK 
and d/ls = 1-6). Therefore the values of the parameters 
Vb, lo and ^ determined by fitting the theory to the ex- 
perimental data shown in Figs. [1] [5] and [3] are consistent 
with Eq. (|16p . Moreover, wc find that the exchange en- 
hanced capacitance F is more than 10 times larger than 
the electrostatic value. This means that an interlayer 
bias V ~ 100 pY can cause a density imbalance of 5%. 
We believe that the density imbalance induced by the 
interlayer voltage plays an important role in the recent 
tunneling experiments performed in the countcrfiow ge- 
ometry (see Section [TV)) . 
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FIG. 4: (a) Experimentally measured I-V characteris- 
tics at d/ls = 1.50 and different temperatures T —: 
25,50,75,100,125,150,175,200,225,250,275 mK. (b) I-V 
characteristics obtained using Eq. (|17|) for Vo = 73 /^V, 
7o = 4.9 pA and different values of a, which are fitted at 
each temperature. The obtained temperature dependence of 
a is shown in Fig. (6] The experimental curves are reproduced 
from the experimental data provided by I. B. Spielman et al. 
and reported in Ref. [3^ . 



B. Temperature dependence 

The measured temperature dependence of the I-V 
characteristic is shown in Fig. |4] (a) . Plotted on a sin- 
gle logarithmic scale, the temperature dependence of the 
peak conductance seems to be consistent with Go (x 
exp(-T/To) as pointed out in Ref. Hi (see Fig. [5]). Such 
a temperature dependence is not expected because it dif- 
fers significantly from the temperature dependence of the 
tunneling current in standard Josephson junction a^^'^^ , 
and it seems naively that the small-bias tunneling con- 
ductance would stay finite even at zero temperaturo^i in 
contrast to theoretical expectationsiSr— . 

We argue that the experimentally observed 
temperature-dependence of the I-V characteristic 
can be described using the parametrization of our the- 
ory. Based on the in-plane magnetic field dependencies 
of the small-signal tunneling conductance measured at 
different temperatures^, we can estimate that ^ depends 
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FIG. 5; Temperature dependencies of (a) the height and (b) 
the width of the conductance peak. Dots show the experimen- 
tal measurements and the lines show the theoretical values 
obtained using Vo ~ 73 fiV, lo = 4.9 pA and the temperature 
dependence of a shown in Fig. [6l The experimental data was 
provided by I. B. Spielman et al. and has been reported in 
Ref . m. 



only weakly on temperature. Moreover, we expect that 
Ps and r do not depend strongly on temperature at the 
experimentally relevant temperatures deep inside the 
coherent phase. This assumption could be tested by 
measuring the temperature dependence of u. Therefore, 
we assume that the main temperature dependence of the 
tunneling current must originate from the temperature 
dependence of r^p. Since the tunneling current according 
to Eq. (|17|) depends on only through the parameter 
q;, we will use a as the fitting parameter. 

Fig. |4] (b) shows the theoretical TV characteristics ob- 
tained using Eq. (|17p . The values of Vq and Iq are as- 
sumed to be temperature independent and the parameter 
a is determined at each temperature as the average of the 
two values obtained by fitting the height and width of the 
conductance peak. The parameter Vq — 72i fiY is deter- 
mined from ^ and u estimated in previous section, and Iq 
is fitted to optimize the agreement between the theoreti- 
cally calculated and experimentally measured height and 
width of the conductance peak at low temperatures^. 
The value of Iq = 4.9 pA obtained by fitting the theory to 




0.05 0.1 0.15 0.2 0.25 



T(K) 

FIG. 6: (color online). Temperature dependence of a. The 
dots show the temperature dependence of a obtained by fit- 
ting the theoretical height and width of the conductance peak 
to the experimental results. The black and red lines are 
guides to the eye. The black line is an exponential curve 
a = c- exp(-ro/r) (c « 15.3 and To « 0.69 K), where c and 
To are fitted using the high-temperature (T > 150 mK) exper- 
imental data. The red line shows a = {T/Tof (To ^ 0.33 K), 
where To is fitted using the low-temperature (T < 125 mK) 
experimental data. Inset shows the temperature dependence 
of Q at small temperatures plotted on a double logarithmic 
scale. 



the experimental data is consistent with the estimate for 
Iq ~ 3.4 — 13.6 pA (corresponding to Asas = 5 — 10 /iK 
and d/lB = 1.5) based on Eq. The theoretically ob- 

tained TV characteristics are in quantitative agreement 
with the experimental ones. In particular, as demon- 
strated in Fig. [5J we can quantitatively obtain the tem- 
perature dependencies of both the height and the width 
of the conductance peak by fitting a single parameter a 
at each temperature. Therefore, the approach outlined 
above is justified and we can concentrate on the temper- 
ature dependence of a, which is shown in Fig.[Sl We can 
separate two regimes of temperatures, where the temper- 
ature dependence of a appears to be different from each 
other. As demonstrated in Fig. [51 at small temperatures 
it is conceivable that a obeys a power law, whereas at 
high-temperatures a shows thermally activated behav- 
ior a cx e""^"/"^. These results seem to be in agreement 
with the expected temperature dependencies of the dif- 
ferent physical processes possibly contributing to vortex 
field fluctuations^^"— The large temperature be- 
havior could be attributed to thermally activated hop- 
ping of merons whereas the power-law dependence could 
result from localized low-energy excitationsi^. Impor- 
tantly, our analysis shows that as T — > our parameter 
a — > indicating that — >■ oo. (At lowest temperatures 
ft/V < fcsT" and it is still decreasing reasonably rapidly 
with decreasing temperature.) 

According to Eq. the conductance diverges in the 
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limit r<p — > cx). Therefore, within our approach a diverg- 
ing conductance at small temperatures is still expected 
by extrapolating the temperature dependence of a, which 
is consistent with the experimentally measured temper- 
ature dependencies of the height and width of the con- 
ductance peak. According to this analysis, the appar- 
ent temperature dependence Go oc exp(— T/Tq) of the 
of the small-signal conductance arises because two differ- 
ent physical processes are contributing to the vortex field 
fluctuations at low and high temperatures, respectively. 

We also point out that can, in principle, be deter- 
mined experimentally using a time-dependent interlayer 
voltage. Similarly to the case of Josephson junctions in 
the presence of strong fluctuations^, we expect that an 
ac field induces structures in the I-V characteristic at 
the resonances eV « nhu! (n = 1,2...). These resonant 
structures resemble the Shapiro steps, but they are dis- 
tinct from each other only if the frequency of the ac field 
satisfies cj > allowing independent determination 

of T^p at different temperatures. Alternatively r^p could 
be determined by measuring the frequency dependence 
of the dynamical conductivity. At the smallest temper- 
atures the required frequencies of the ac field arc in the 
GHz frequency range. 



IV. QUALITATIVE CHANGES CAUSED BY 
LARGER TUNNELING AMPLITUDE 

We next study the effects caused by a larger tun- 
neling amplitude. In a series of recent tunneling 
experiments^^— samples with otherwise similar param- 
eters but significantly larger tunneling amplitudes were 
used. The tunneling amplitude can again be determined 
either by using the tunneling I-V characteristics in the 
absence of magnetic field or by solving the Schrodinger 
and Poisson equations self-consistently. The obtained 
values of As as vary between 100 - 150 ^Kiii^. This 
good agreement between the values obtained using two 
different methods suggests that the estimates for A^yis 
arc reliable. In comparison with the experiments dis- 
cussed in the previous section, the tunneling amplitude 
of the samples in these experiments is therefore approx- 
imately an order of magnitude larger. Based on the 
scaling law Iq oc Ag^g, Eq. ([TO]) , we expect that the 
tunneling currents are approximately two orders of mag- 
nitude larger Iq ^ 0.5—1 nA^. This kind of scaling 
is in agreement with the experimental observations^^—. 
Moreover, it was observed in the experiments^ that the 
tunneling current is still proportional to the area of the 
sample. These observations indicate that essentially the 
same fluctuations-dominated physics determines the tun- 
neling I-V characteristics also in these experiments. On 
the other hand, there are several important qualitative 
differences, which are discussed below. These new effects 
arise because the tunneling resistance is no longer large 
compared to the other resistances in the system. 



A. Tunneling geometry 

Because Iq is now two orders of magnitude larger than 
in the earlier experiments, we estimate that the tunnel 
resistance at low temperatures is comparable to the con- 
tact resistances and intralayer resistances Rxx and Rxy 
Although the interlayer voltage V{r) is no longer con- 
stant, we can still expect that it varies reasonably slowly 
in space. This means that we can describe the trans- 
port in this system by assuming that the local tunnel 
current density is I{V{f))/L'^, where V{r) is the local 
interlayer voltage and I{V) is the current-voltage char- 
acteristic given by Eq. ([TT]). The interlayer voltage V{'r) 
can now be calculated by assuming that the intralayer 
transport and counterflow currents are described by re- 
sistivity tensors obtained from the experiments^i^i^"— . 
An effective transport theory under these assumptions 
can be formulated using the current continuity equation 
and appropriate boundary conditions, which arc deter- 
mined by the experimental geometry. 

Here we do not attempt to explain quantitatively all 
the details of the experiments but instead we just want to 
qualitatively explain the main features observed in the I- 
V characteristics. The main effect of the inhomogeneities 
in is that all the sharp features in the I-V character- 
istics become rounded and dependent on the locations of 
the voltage probes. In particular, it becomes difficult to 
determine the tunneling conductance at small voltages, 
because the tunneling resistance is small compared to the 
other resistances in the system and the observed values 
of the interlayer voltage can depend on the locations of 
voltage probes^. Nevertheless, the temperature depen- 
dence of the peak tunneling curren t ^"^'^^ is similar to that 
observed in the earlier experiments (Section IIIIBp and 
therefore we expect that the perturbative approach is at 
least approximately valid even at the lowest experimen- 
tal temperatures. By keeping in mind the type of small 
modifications discussed above, we can mainly ignore the 
inhomogeneities in the analysis of the remaining effects 
caused by larger tunneling amplitude. 

FigS-EUa), (c) and (e) show the results of experimental 
measurements at d/ls ~ 1-42 and T < 20 mK, reported 
in Ref. H^. The I-V characteristic shown in Fig. ^c) is 
qualitatively very similar to the I-V characteristics ob- 
tained for samples with smaller tunneling amplitude (see 
Fig. [3]) . The main difference to the earlier experiments 
is that there now exists jumps at particular values of the 
voltage. If we assume that the inhomogeneities are not 
important, the experimental geometry can be described 
with an effective circuit shown in the inset in Fig. [7] (b). 
This gives us a circuit equation 

V2t -RI^ V{I), (18) 

where the V{I) characteristic of the bilayer tunnel junc- 
tion is obtained from Eq. (|17p by solving the voltage V 
as a function of current /. The solutions of Eq. (fT5|) are 
obtained by finding the intersections of V{I) characteris- 
tic and the lines V2t — RI ■ We assume that a = 0.01 and 
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FIG. 7: (color online). (a),(c) Experimentally measured de- 
pendencies of the interlayer voltage V and the tunneling cur- 
rent / on the applied voltage Vit- (e) Experimentally mea- 
sured tunneling current / as a function of interlayer voltage V . 
The experimental measurements were made at d/Zs = 1.42 
and T < 20 mK. (b), (d), (f) Corresponding theoretically 
calculated tunneling characteristics obtained using Eqs. (|18|) 
and (|17|l . The theoretical parameters used in the calculation 
are a = 0.01, Vb = 100 pN and Jo = 0.56 uA. The inset in 
figure (b) shows the effective circuit. Resistance R is caused 
by the circuit and contact resistances. There exists a range 
of voltages V2t, where Eq. (|18p has several solutions. The 
additional solutions are indicated by the thin dashed lines in 
figures (b), (d) and (f). The experimental curves are repro- 
duced from the experimental dataprovided by L. Tiemann et 
al. and reported in Refs. 13 and |45| . 



Vb = 100 in agreement with the values found in Sec- 
tion lllll Because the effective resistance R is significantly 
larger than the tunneling resistance at small voltages, the 
resistance R can be determined from Fig. ^c) by calcu- 
lating the slope of the I{V2t) characteristic at small ap- 
plied voltages V2t- We obtain R « 150 kfi. Moreover, by 
fitting the theoretically calculated maximum of the tun- 
neling current to the experimentally observed maximum 
current, we obtain Iq = 0.56 nA. This value is in good 
agreement with an estimate /q ~ 0.65 nA obtained using 
Eq. (Ull) for d = 28.6 nm, d/h = 1.4, L^Ly = 880 • 80 
firn^, A5^5 — 100 /iK and ^ = 92 nm (corresponding to 



Vq = 100 /iK). Here we have assumed that u, ps and F 
are the same as in the Caltech experiments. The results 
of the theoretical calculations for these values of a, R 
and /q are shown in Figs. [7] (b), (d) and (f). They are 
in good agreement with the experimental measurements. 
Moreover, we find that for particular values of V2t, there 
exists several solutions of Eq. (|18p indicating bistability. 
We expect that the jumps in the measured I-V character- 
istic [see Fig. [7] (e)] originate from a hysterctic switching 
caused by this bistability. 



B. Counterflow geometry 

Recently, tunneling experiments were performed also 
in the counterflow geometry [Fig. [5] (a)] using samples 
with similar As^s ~ 100 — 150 /iK4^ In this case, dc 
input current is generated by applying a voltage to a 
large resistor in series with the system, and the current 
which passes through the top layer is redirected via a 
loop resistor into the bottom layer. If we assume that 
inhomogencities are not playing an important role, we 
can construct an equivalent circuit shown in Fig. [Sl^b). 
The circuit equation is now 



'Total 



= /Tunnel + /Loop = I{V) + V / R, 



(19) 



where liV) is the tunneling I-V characteristic of the bi- 
laycr and R is an effective resistance, which takes into 
account the loop resistance Rl and the contact resis- 
tances. 

Two intriguing observations were made in the experi- 
ment reported in Ref . |4^. It was found that the tunneling 
current was large as long as the injected drive current was 
smaller than a certain critical current Ic at which point 
a transition into a new regime with negligible tunneling 
current was observed. For small (I/Ib the observed tran- 
sition is extremely sharp and becomes smoother for larger 
(I/Ib as shown in Figs. [5] (c) and (d), respectively. Sec- 
ondly, it was found that although the transition between 
the two tunneling regimes had only a minor impact on the 
state as observed in the magnctotransport, the longitu- 
dinal resistance increased rapidly for currents exceeding 
the critical current Ic- Our theory provides a simple and 
natural explanation for both experimental observations. 

In order to explain these observations we assume that 
a = 0.01 and Vq = 100 /iV similar to the values used in 
previous sections. We determine R = 180 kfl and Iq = 
0.16 nA {d/lB = 1.68) and Iq = 0.48 nA {d/lB = 1.37) 
such that we visually obtain good agreement between 
theoretical and experimental dependencies of loop and 
tunneling current on the total current. The values of Iq 
are in agreement with the measured tunneling currents 
at different d/Zs in the tunneling geometry^i^. The ef- 
fective resistance R is significantly larger than the loop 
resistance Rl = 10 kfi used in the experiment^, and 
therefore we expect that small variations in i?^ would 
not affect the tunneling and loop currents significantly. 
Our calculations show that depending on the value of Iq 
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FIG. 8: (color online) . (a) Counterflow geometry, (b) Equiva- 
lent circuit assuming homogeneous interlayer voltage. Here R 
is an effective resistance, which takes into account the loop re- 
sistance and contact resistances, (c), (d) Tunneling and loop 
currents as a function of total current. The experimental re- 
sults shown in figures (c) and (d) are obtained for d/lB = 1.37 
and d/ls = 1.68, respectively. In theoretical calculations we 
assume that R, a and Vo are independent of d/Zs. The param- 
eter 7o is fitted to be Jo = 0.48 nA and Jo = 0.16 nA for the 
different values of d/ls = 1-37 and d/ls = 1.68, respectively, 
(e), (f) Interlayer voltages corresponding to the parameters 
used in figures (c) and (d), respectively. For lo — 0.48 nA 
there exist an interval of Ixotai, where the Eq. 1)19^ has sev- 
eral solutions. The additional solutions are shown by the thin 
dashed lines in figures (c) and (e). The experimental curves 
are reproduced from the experimental data provided by Y. 
Yoon et al. and reported in Ref. [i^ . 



there exist two qualitatively different behaviors similarly 
to the ones observed in experiments. By inspecting the 
tunneling I-V characteristics (Figs. [Hand [7]), we see that 
after the peak tunneling current has been reached the 
tunneling current decreases quickly with increasing volt- 
age. If /o/Vq is large enough this decrease is more rapid 
than the increase of the loop current through the resis- 
tor. In this case, there exist several solutions of Eq. 
around the critical current as demonstrated with the help 
of dashed hues in Figs.|S](c), (e). We expect that this re- 
sults in bistability, which can be seen as abrupt changes 
in the tunneling and loop currents in the experiments at 
small d/ls- At large (I/Ib the ratio Iq/Vq is expected to 
be smaller and therefore the change in the tunneling and 
loop currents becomes smoother near the critical current. 
Our theoretical picture explains naturally why the tran- 
sition has only a minor impact on the ut = ^ state as 
observed in magnetotranspori4^. The reason is simply 
that the transition occurs due to the tunneling I-V char- 
acteristic, which is an intrinsic property of the vt ~ ^ 
state, and therefore in the first approximation it should 
have no effect at all on the 1/^ = 1 state. Moreover, ac- 
cording to our theory it is natural that the critical current 
observed in this counterflow geometry is almost exactly 
the same as the peak tunneling current observed in the 
tunneling geometry. 

As already pointed out, the experiments also indicate 
that the dissipation increases rapidly at the currents ex- 
ceeding the critical current. It has been suggested^ 
that the increase of dissipation could be caused either 
by meronsiiri£iiSi2S or so-called Josephson vortices^i. 
However, the Josephson vortices could only exist if the 
Josephson length Aj = AttP^Ps/ l^sAS was smaller or at 
least comparable to the coherence length ^. This is not 
the case in the present experiments where we estimate 
Aj ^ 4 /im. Moreover, the existence of Josephson vor- 
tices would indicate significantly larger critical currents^! 
than observed in the experiments. Therefore we expect 
that merons are responsible for the increased dissipation. 

In our theoretical picture, the increase of dissipation 
should be due to a change in the dynamics of merons, 
which can be due to either an abrupt change in interlayer 
voltage or due to large counterflow current. If the effect 
was caused by the counterflow current it seems difficult 
to understand why the critical current is almost exactly 
equal to the maximum current observed in the tunneling 
geometry. Moreover, by assuming that the counterflow 
current induces a Magnus force on meronsi^ and that 
merons are confined into potential wells corresponding 
to the measured thermal activation gap^, we estimate 
that the critical counterflow current necessary to depin 
the merons is significantly larger than typical experimen- 
tal currents. On the other hand, we find that for small 
d/ls the interlayer voltage significantly increases at the 
critical value of total current. Moreover, for both val- 
ues oi d/ls, the interlayer voltage increases rapidly with 
increasing total current for currents exceeding the criti- 
cal current [Figs. [5] (e), (f)], resembling the experimen- 
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tally observed dependence of the dissipation on the total 
current^. Therefore, we believe that the increasing dissi- 
pation is tightly connected with the increasing interlayer 
voltage. By using our estimate for the exchange enhanced 
capacitance F, we find that an interlayer voltage of 100 
/iV can result in a density imbalance of several percents. 
It has been experimentally observed^ and theoretically 
explainedi^ that such a density imbalance can signifi- 
cantly affect the longitudinal resistance. Furthermore, it 
is possible to experimentally test this hypothesis, because 
the density imbalance would result in different longitudi- 
nal resistances in the upper and lower layersi^i^. 



V. RELATION TO PREVIOUS THEORETICAL 
APPROACHES 

The theoretical explanation of charge transport in bi- 
laycr quantum Hall systems has been a long-standing the- 
oretical problem. Perhaps most surprisingly, instead of a 
true Josephson effect, a Josephson-like conductance peak 
has been observed experimentally. Understanding the 
finite height and width of this small bias conductance 
peak and the calculation of critical tunneling currents 
has attracted a significant amount of interest. Theoret- 
ical approaches differ in whether they take a clean or a 
disordered system as a starting point of their analysis. 

In the clean limit, the critical current is determined as 
the maximum current for which the pseudospin field can 
be static, and the dynamics of charge transfer between 
leads and bulk by quasi-particlcs plays an important role. 
Conductance peak heights and maximum tunnel currents 
predicted in the clean limit are orders of magnitude larger 
than the values observed in experiments?^ and also dif- 
fer from experiments with respect to the parametric de- 
pendencies of conductance peak and maximum tunnel 
current on sample area, bare tunneling amplitude, and 
in-plane magnetic field. Most notably, the critical tunnel 
current depends linearly or sublinearly on the tunneling 
amplitude AsaSi and is independent of system size for 
systems larger than the Josephson length Aj. In contrast, 
experiments seem to be consistent with a critical tunnel 
current proportional to both the area of the sample and 
to the square of the tunneling amplitude. 

The influence of quenched disorder can be described 
by the introduction of a vortex field, which causes an ex- 
ponential decay of pseudospin correlations with a char- 
acteristic spatial correlation length ^. Under the as- 
sumption of a static vortex field, local pseudospin cor- 
relations are assumed to have an infinite range in time. 
One-dimensional numerical simulations indicate that the 
static vortex field results in an intriguing self-organized 
Bean critical state, where the current injected at the 
boundaries can affect the condensate dynamics through- 
out the whole sample.— By generalizing this finding to 
two-dimensional system, the authors in Refs. and [2^ 
find that if the sample size is sufficiently large, the crit- 
ical tunneling current is proportional to the area of the 



sample. This statement can be justified with the help 
of an energy optimization argument in the presence of 
random static vortex field^ii^, and it turns out that the 
sample size should be compared to a length scale Xj/S,- 
In the experiments reported in Ref. 13 the conditions 
Lx,Ly > A J are approximately satisfied and the criti- 
cal tunneling current proportional to the area of the sam- 
ple has been observed. On the other hand, the approach 
based on the assumption of a static vortex field does not 
easily explain why a similar scaling with the area has 
been observed^ also in samples, where the sample size is 
not large enough in comparison to the Josephson length. 
Other challenges for this type of approach are the expla- 
nations of the experimentally observed temperature and 
parallel magnetic field dependencies of the tunneling cur- 
rent, and the quantitative description of the height and 
width of the conductance peak. Nevertheless, one would 
expect that the vortex field should be static in the limit 
of zero temperature, and that the predictions of Refs. [U 
and [2^ should describe the limit of extremely low tem- 
peratures, which probably has not yet been reached ex- 
perimentally. 

Effects of a dynamical vortex field were analyzed in 
Refs. [lOl and [m. In the present manuscript, we follow the 
approach outlined in Ref. [Tol . introduce a phenomenolog- 
ical vortex field, and calculate the dependence of exper- 
imental obscrvablcs on the correlation length ^ and the 
correlation time of the vortex field. We argue that all 
the important temperature dependence of experimental 
observables originates in the temperature dependence of 
the dimensionless decoherence rate a{T) = £^/{ut^{T)), 
and we determine a(T) by fitting experimental curves. 
The dynamics of merons was explored in Ref. [13 by 
means of mapping the quantum Hall bilayer to a classical 
two-dimensional XY model with a symmetry-breaking 
field and with disorder. It appears that disorder nucleates 
strings of overturned spins, which connect vortices and 
antivortices at their ends. At low temperatures, this state 
is characterized by anomalously large fluctuations of the 
vortex field, which are reminiscent of glassy dynamics. 
We believe that the power law like, weak temperature 
dependence of a which we extract from experiment may 
be consistent with the results of Ref. [l^. 



VI. DISCUSSION AND CONCLUSIONS 

In this paper we have studied how the Josephson- 
like tunneling depends on the area and tunneling am- 
plitude of the samples, applied parallel magnetic field 
and temperature. We have compared different theoretical 
approaches describing the Josephson-like tunneling and 
found that all the tunneling experiments are in agreement 
with a theory which treats fiuctuations due to meron ex- 
citations phenomenologically and takes tunneling into ac- 
count pcrturbativcly. Previously, Josephson-like tunnel- 
ing in samples with smaller— i^^i^^—i^"— and larger—"— 
tunneling amplitudes have often been discussed in differ- 
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ent terms, because two orders of magnitude larger tun- 
neling currents are found in latter experiments. Never- 
theless, our analysis shows that all these experiments can 
be discussed within a common framework, and the ob- 
served magnitudes of the tunneling currents are in good 
agreement with the theory. We have also shown that 
the somewhat surprising temperature and in-plane mag- 
netic field dependencies of tunneling current can be ex- 
plained quantitatively by describing the effects caused by 
the merons with the help of a phenomenological vortex 
field. In this paper, we have concentrated on the de- 
scription of the tunneling deep inside the coherent phase 
allowing us to make several simplifying approximations 
in Section |lll However, we note that in principle the 
equations (fTOj) . (fTTj) . (fT2|) and dH) could also be used to 
describe the Josephson-likc tunneling closer to the phase- 
boundary separating the coherent and incoherent phases. 
This would allow the determination of the theoretical pa- 
rameters of the model function of different exper- 
imentally controllable quantities such as d/ls, density 
imbalance and temperature. 

The surprisingly good agreement found between the 
theory and experiments naturally introduces several new 
theoretical questions. The temperature dependence of 
the tunneling currents was explained by fitting the pa- 
rameter , which describes the dynamics of the vortex 
field. However, the details related to the microscopic ori- 
gin of the vortex field dynamics are not well-understood. 
The exponential temperature dependence at large tem- 
peratures most likely results from the thermally activated 
hopping of the merons, but the mechanism giving rise 
to the power law dependence at low temperatures re- 
quires additional theoretical analysis. Other interesting 
questions are related to the validity of the main approx- 
imation in our approach, which is the treatment of the 
tunneling as a perturbation. In general, the validity of 
this approximation depends on the nature of vortex-field 
fluctuations and dissipation mechanisms. We do not try 
to specify the dissipation processes here. Instead we ex- 
pect on general grounds that the pcrturbative approach 
is controlled by a condition 



condition 



V h' 



(20) 



where defines the size of a domain which is large 
enough such that its dynamics can be considered to be 
independent of the other domains. The idea behind this 
inequality is that the left side defines the conductance of 
an independent domain and we interpret it as a transmis- 
sion probability multiplying /h. If the condition ((20)) 
is satisfied, the transmission probability is much smaller 
than one and the higher order tunneling processes are 
suppressed by the powers of the transmission probability. 
There are two relevant length scales, ^ and ur^p, in the 
problem, which could determine the size of the indepen- 
dent domain Li. Because at low temperatures wr^ ^ ^, 
the assumption L| sa {ut^Y certainly yields the sufficient 



ire 



-Do/2 



< 1 



a"- 



(21) 



for the applicability of perturbation theory. In Appendix 
IB| we obtain another condition (jBlSP for the validity 
of the pcrturbative approach in terms of the dissipa- 
tion strength 7. For dimensional reasons, we expect that 
(X l/r^ so that the only possible difference between 
conditions (|2T|) and (jB15|) is the reasonably large pref- 
actor ps/{h/Tip) appearing in condition pip. This pref- 
actor essentially originates from our assumption for the 
area of the independent domains. The result of the 
classical calculation in Appendix [B] can be interpreted as 
suggesting that the area for independent domains should 
be chosen as Lif « wr^^, giving a less restrictive criterion 
for the applicability of perturbation theory. 



Because the assumption used to obtain Eq. ([2T|) likely 
overestimates Li, we consider the condition (|2ip as a suf- 
ficient condition for the validity of the pcrturbative treat- 
ment of the tunneling, and we are confident that if it is 
satisfied, the pcrturbative approach is well justified. We 
find that the condition (PT|) is always satisfied in the ex- 
periments with smaller tunneling amplitude. In the case 
of larger tunneling amplitude, assuming that the tem- 
perature dependence of a does not change, this condi- 
tion could break down (depending on the exact values of 
different theoretical parameters) at the smallest experi- 
mental temperatures. The experimental observation^^— 
that there is no clear transition to a different operation 
regime as a function of temperature could indicate that 
the necessary condition for the pcrturbative treatment 
of tunneling is less restrictive than the condition (pi]) or 
that the larger tunneling currents resulted in an increase 
of the dissipation and the fluctuations. Experiments with 
even larger Ag^s could reveal this type of transition. 

In this paper, we have discussed the Josephson-like 
tunneling at small voltages, where the theory and ex- 
periments are in good quantitative agreement with each 
other. On the other hand, it is clear (see Figs. [1] and 
S]) that our theoretical predictions differ from the ex- 
perimental results at larger voltages. This discrepancy 
could be caused by the higher order terms in the Hamil- 
tonian [Eq. a slightly incorrect form of the vortex 
field correlation function [Eq. ([9])], the effects caused by 
the dissipation (sec Appendix[B]) or an incoherent tunnel- 
ing occurring inside compressible puddles of the electron 
liquid^^. The last possibility is supported by the exper- 
imental observation that at significantly larger voltages 
the tunneling I-V characteristic starts to look more and 
more similar to the corresponding I-V characteristic in 
the incoherent phase^i^. Theoretical treatment of the 
tunneling at large voltages goes beyond the scope of the 
present paper. 
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Appendix A: Calculation of the I-V characteristics 

Similarly as in Josephson junctions^ we can use 
Fermi's golden rule to calculate a forward tunneling rate 



^ieVt/h^iQB{xi-X2) 



X (^giv{ri,t) ^ifmiri.t) ^-iip{r2,0) ^-iipm{r2,0)'^ (Al) 



Due to the symmetry of the problem the backward tun- 
neling rate is related to the forward tunneling rate as 

^iV,B)=f{-V,-B). (A2) 
By assuming that if and ipm commute, we get 



the commutation relations are [7rjr,((5j,] = 
the Hamiltonian defined by Eq. ([S]) can be rewritten as 



hSj: and 



H = / (frn 



2AI ^ 



Mujl 

k 



2M 



(A8) 



where M = Y?V je^ and a;^: is the dispersion relation for 
the collective Goldstone mode given by Eqs. (jG]) and 

The creation and annihilation operators and at are 
defined as 



4 



2h 



(A9) 



It is easy to show that they satisfy commutation relations 

[oj:, at,] = '^fcfco ['^L'^l'l ^ ^' i^k'^ic'] ~ *-* ^^'^ allow to 
rewrite the Hamiltonian as 

H = + hurl aiar + 1/2 



(AlO) 



We can now express ip{f,t) as a Fourier series 

(AU) 

where the summation is taken so that each pair ±fc is 
counted only once. Thus the correlation function can be 
written as 



/ = e[f{V,B)-fiV,B)] 
2-Ke\^L^ 



■[S{QB,eV)-S{- 



~eV% (A3) 



where 

S{QB,eV) 



1 

2TTh 



dt / d^r e'"^*/'^e*'5«^ 

> J 

X (e*^("^'*)e-^^(°'°))G„(r, t). (A4) 



and 



G^{r,t) = (e'v™(--*)e-*V'"(o^o)). (A5) 



;[¥>£(t)e"'-'- + tp_g(t)e- 



(A12) 



where the product is taken similarly over the pairs ±fc to 
avoid double counting. 

By using the Baker-Hausdorff relation, we can rewrite 
Eq. (|XT2)) as 



g-iC(f,t)/2/ f [.<5o(t)-¥'o(0)]\ 



(A13) 



By defining the Fourier components of the operator 
fields (p and tt as 



(A6) 
(A7) 



where 
C(f,t) 



I 

l2 



(A14) 
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By using the Wick theorem for equilibrium correlation 
functions, it is easy to show that an operator ip, which is 
a linear combination of creation and annihilation opera- 
tors, satisfies 

^(e'"^) = -a(V''>(e*"'^), (A15) 
da 

and therefore {e^"^') — e^^^'"^'^ . 

By applying this result to Eq. (|A13p . we can rewrite 
the correlation function as 



Jv{r,t) -i<f {0,0) \ 



~tC{r,t)/2~D(r,t)/2 



(A16) 



where 
D{f,t) 



1 



/3ML2 L2 



The functions C and D can now be calculated by ex- 
pressing Lpj^ in terms of creation and annihilation opera- 
tors. After a straightforward calculation, we obtain 

Cir,t) = ^|^+E^-K^^)-^(fc-^| 
1 fe^t uh . , , /r ^xl 

fe#0 

= -^^^sin(t^ji)cos(fc-r). (AI8) 

Note that in each step we have changed the summation 
convention, so that finally the summation is over all k 
vectors. In the final form the term corresponding to k ~ Q 
should be understood as a limit fc for the given 
expression. 

It is easy to show that C is independent on the direc- 
tion in the x,y-plane C(f, t) = C(r, t) and that C(r, —t) ~ 
-C(r,t). Moreover, for t > 



-e{ut~r) 



1 



(A19) 



t2 - 



The last expression is easiest to prove by calculating the 
Fourier expansion for the right side of Eq. (IA19[) . 
Similarly, we get 



D{r,t) = Y. 



hu 



L^Psk 

k 



[1 — cos(a;rri) cos(fc • r)] coth{/3hidr/2), 



(A20) 

where the fc = term should be calculated by setting 
cos(fc-r) — > 1 and taking the limit /c — > for the rest of the 
expression. It is again easy to show that D{f, t) = I?(r, t) 
and D{r,-t) = D{r,t). 

By using Eqs. (|X3)) . (|X4)) and (|A16p and the properties 
C{f,t) = C{r,t), D{r,t) ^ D{r,t), C{r,^t) = -C{r,t), 
D{r,—t) ~ D{r,t) and G„i{r,—t) = Gm{r,t), we get 
equations (HU]), (HI]) and P^ . 



Appendix B: Derivation of the I-V characteristics 
using Sine-Gordon equation 

The Lagrangian density corresponding to Hamiltonian 
density ([5]) is 



C 



2e2 



if'' - ^{VipY + 2XC0S 



^P + ^Pm +QbX 



eVt 
(Bl) 

The equation of motion for ip can be written as a Sine- 
Gordon equation 



1 



V </3 ^9* V - 1"^ 



1 



P + Pm + Qbx + eVt/h 



(B2) 

where we have introduced an operator 7, which describes 
the dissipative terms in our model. The actual mecha- 
nism of the dissipation is not well- understood, but several 
different dissipative terms have been proposedi^i^'2^. 
Here we keep the dissipative terms unspecified for as long 
as possible. 

The idea is to formulate the perturbative tunneling 
theory so that we would only need to assume that the 
variations of p are small within certain reasonably large 
space and time scales. More concretely, we divide the 
space-time into domains (labeled with index i) of size 
Ti (timc-intcrval) and (spatial area). We denote the 
average value of p inside the domains as p^, and as- 
sume that the variations within the domain are small, 
5pi{ri,ti) = p>{fi,ti) — p^'^l. The variations Spi inside 
each domain i are not entirely independent on the other 
domains because the domains can be dynamically cou- 
pled with each other. Here we assume that the domains 
are large enough that the boundary conditions are not 
important. Although the boundary conditions strongly 
affect the state of the whole system in the case of static 
vortex field^l, we believe that our approximation can be 
justified in the typical experimental situation where the 
fiuctuations of the vortex field are important. 

Inside each space-time domain we need to solve an 
equation 



^2 r d?Spi 1 . 

V dp^i 7d(/5j = Tj^^^ 

u \j 



eVt 

Pi+iPm+QBX+^- 



(B3) 

The time average of tunneling current liV) can then be 
calculated from equation 



dti 



x(sin [5pi{f,,ti) + Pi + Pm{n,U) + QuXi + eVti/fi\). 

(B4) 

Here — >■ 00 is the total time, (...) denotes the en- 
semble average over different realizations of the vortex 
field and the integrations are over the space-time do- 
mains. 
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It is easy to show that a particular solution of Eq. (jB3[) 
can be written as 



1 



(fir' dt' Gin - r' ,U ~ t') 



X sm 



Tp, + Lprair', t') + Qbx' + eVt'/h 



(B5) 



where 



L2 27r 



(B6) 



and 



'xG{fi — r' ,ti ~ t') sin 



(fir dt' 

Tp^ + ipmifi, t') + Qbx' + eVt'/h 

(B7) 



is a constant guaranteeing that ip^ is the average value of 
(p within the domain i i.e. 



dU / d^ri5ipi{ri,ti) = 0. 



(B8) 



We have also assumed that the dissipation operator act- 
ing on the plane waves satisfied the following equation 

-gifc.r-g-ic^t ^ ky^-^e-'^^K (B9) 

We assume that the dissipation is isotropic in the (x,y)- 
plane. 7(0;, A:) = j(u),k). Moreover, we require that 
G(— cj, — fc) = G*{uj, k), which gives 7(— w, k) = —j{uj, k). 

By assuming that the instantaneous value of the vortex 
field pm {r, t) and the average quantities Gi and Tp^ are 
statistically independent, we obtain using Eqs. (jB4|) and 
(IBSP that the tunneling current liV) can be written as 



1 e As AS 1 



2 /i 47rZ| A2 



Im 



,i[¥'m('",t')+Qi3a;'+eVt7?i-v™(ri,t)-Qi5a;i-eyti/;i] 



(BIO) 

Here we have taken into account that the vortex field 
tends to stay constant over some time and space inter- 
val and therefore the differences of the vortex field val- 
ues ipmir',t') — ipm(fi,t) are correlated. On the other 



hand, the actual value of the vortex field ipmir, t) can be 
arbitrary, and therefore only the terms proportional to 
^±i[^Pm{r ,t )-ip„,(ri,t)] contribute to the average cur- 
rent. Physically this means that the fluctuations of the 
vortex fleld destroy the supercurrent. 

By assuming also that correlation function for the vor- 
tex field is given by Eq. we obtain [R = r/^, q = fc^, 

n = u;^/u, y{n,q) = ^{nu/^,q/o] 



(!) = 



I'. 



dq / dii 

Jo 



y{n,q)e 



2q 



h'in,q)e? 

a a 

+ {-n + v/Vo)^ ^ a^ + {n + v/Vofi 



X j dRRe-^Jo{QB^R)MqR)^, (Bll) 

where Vb and a are given by Eqs. ([J) and ([3]), respec- 
tively. The current scale is now determined by 



which differs from the expression (IT6l) by a factor e -^o/^. 
If we further assume that ^i'{Vl,q)S^^ <C 1, we can use 



lim — - 



1 



to obtain 

(/> ^ rJdq 



[% + yo) + ^(y-2/o)] 

(B13) 



cfi + {V/Vo - qfi o? + (y/^o + q? 
xJdRRe-"JniQBCR)M<lR)- (B14) 

The only difference to the expression ([T5|) is that the 
current scale Iq determined by Eq. (jl6p is now replaced 
by given by Eq. ((BT2|) . We see from Eq. (|BTT|) that 
in the general case also the dependence of the dissipation 
on the frequency and momentum of the pseudospin waves 
can affect the I-V characteristic. 

In order to determine the limits of validity of the 
perturbation theory we need to calculate {[ip(r,t) — 
^(0,0)]2) - M0,0)]2) = 2(V(r,i)) - 

2{6p{r,t)Sip{0,0)). Using similar methods as above, we 
obtain that the perturbation theory is controlled by a 
condition 



r 1 



(B15) 



where for simplicity we have assumed that the dissipation 
7 does not depend on frequency and momentum. 
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There are some sample to sample variations also in other 
parameters and therefore it is not meaningful to try to 
make the prediction for 7o more accurate. However, we 
point out that it would be very useful to study experi- 



mentally the scaling of Jo as a function of As as by using 
samples which have otherwise as identical parameters as 
possible. 



